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Time  Durations  of  Rain  Rates 
Exceeding  Specified  Thresholds 


I.  INTRODUCTION 

The  attenuation  of  electromagnetic  energy  caused  by  rainfall  has  posed  prob¬ 
lems  in  the  operational  efficiencies  of  many  weapon -guidence,  surveillance,  and 
communication  systems  and  continues  to  be  a  concern  in  those  currently  in  design 
or  development.  Depending  upon  the  intended  use,  each  such  system  operates  on 
a  specific  wavelength  or  series  of  frequencies  that  may  show  different  degrees  of 
attenuation  with  various  intensities  of  rainfall.  Different  rain  conditions  may  re¬ 
quire  the  choosing  between  systems  and/or  frequencies  to  insure  effective  opera¬ 
tion.  The  ability  to  predict  the  performance  of  a  particular  system  in  any  antici¬ 
pated  rain  situation  would  thus  be  advantageous. 

A  performance  evaluation  of  any  system  requires  a  two-step  effort  with  the 
first  being  the  definitive  knowledge  of  the  attenuation  that  is  experienced  from 
using  particular  frequencies  in  respect  to  the  absolute  rain  rates.  This  informa¬ 
tion  may  then  be  associated  with  real-time  measurements,  forecast  projections, 
or  past  climatological  records  of  operational  areas  to  predict  periods  and  degrees 
of  attenuation. 


(Received  for  Publication  22  May  1985) 


The  attenuation  of  a  particular  frequency  can  be  determined  by  the  direct 
measurement  of  the  electromagnetic  energy  that  is  propagated  through  rain  of 
known  rates  over  prescribed  path  lengths.  Since  rain  intensities  may  vary  con¬ 
siderably  within  short  time  periods,  detailed  records  of  rain  rates  along  the  line- 
of-sight  are  essential  for  accurate  attenuation  determinations.  This  is  not  an 
insunn ountable  problem  however,  now  that  instruments  are  available  for  the  high- 
resolution  measurements  of  rain  rates.  One  such  device  was  designed  and  con¬ 
structed  at  the  Air  Force  Geophysics  Laboratory  and  a  detailed  description  of  the 
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instrument  and  the  techniques  used  in  the  data  reduction  have  been  reported.  ' 

The  more  immediate  problem  is  the  prediction  of  rain  rates  on  which  to  base 
attenuation  assessments.  The  most  idealized  scenario  would  be  a  perfect  forecast 
of  the  time -resolved,  rate  structure  of  an  impending  rain  situation.  In  this  case, 
the  times  where  attenuation  becomes  a  problem  could  easily  be  defined  by  noting 
when  the  rain  rate  reaches  or  exceeds  a  critical  intensity.  This,  of  course,  is 
beyond  the  ability  of  any  present  forecasting  technique,  whether  it  be  by  humans 
or  computers.  Realistically,  the  best  that  can  be  provided  now  are  estimations  of 
total  rainfall,  time  duration,  and  some  possible  rate  intensities. 

Climatological  records  are  available  for  most  operational  areas  and  they  pro¬ 
vide  past  information  on  the  number  of  rain  situations,  total  and  mean  amounts, 
and  rain  rates,  although  the  listed  rate  measurements  would,  more  than  likely,  be 
averaged  values  over  long  time  periods  of  5  min  to  1  hr  or  more.  These  statis¬ 
tics  may  be  used  to  provide  gross  estimates  of  attenuation  but  until  a  reasonably 
accurate,  rain-rate  prediction  method  is  devised,  attenuation  assessments  will 
remain  inadequate. 

This  report  describes  our  investigation  into  rain-rate  variability.  Our  data 
base  consisted  of  31  days  of  high-resolution,  rain-rate  measurements  taken  at 
Hanscom  Air  Force  Base  (HAFB)  during  1982  and  1983.  The  study  was  approached 
from  the  viewpoint  of  the  operational  commander  who  has  to  predict  the  efficiency 
of  his  system  from  forecast  information  and  that  of  the  engineer  who  has  to  design 
the  most  efficient  system  for  use  in  some  specified  area  whose  weather  history 
can  be  found  in  the  climatological  records.  Both  concerns  can  be  reduced  to  one 
common  question,  "How  often  and  for  how  long  will  a  particular  system  suffer  at¬ 
tenuation  while  operating  in  a  rain  situation  that  can  only  be  described  by  forecast 
or  climatology?". 
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2.  DATA  DESCRIPTION 


The  31  days  of  rain-rate  measurements  that  were  taken  at  HAFB  during  the 
summer  and  fall  seasons  of  1982  and  1983  consisted  of  periods  of  rain  that  varied 
from  a  few  minutes  to  ~  5  hr.  The  maximum  recorded  rate  was  84  mm  hr  *. 

It  was  decided  that  the  data  that  were  to  be  used  in  the  analysis  would  only  in¬ 
clude  those  periods  when  it  was  actually  raining  measureable  amounts.  Thus,  a 
trace  of  rain  having  a  rate  of  0.  127  mm  hr  *  (0.  005  in.  hr  *)  was  designated  as  a 
lower  limit.  After  applying  this  limitation,  the  data  were  divided  into  two  cate¬ 
gories  with  (1)  numerous  instances  of  rain  having  durations  less  than  1  hr  and 
(2)  16  instances  of  rainfall  lasting  longer  than  1  hr. 

In  the  first  category,  no  rain  instance  was  considered  for  time  periods  of 
<  5  min.  These  cases  had  fewer  points  for  statistical  evaluation  and  were  invari¬ 
ably  those  of  extremely  light  rain  and  it  was  found  that  those  particular  conditions 
were  amply  represented  by  numerous  cases  of  >  5-min  duration.  In  the  second 
category,  the  rain  data  were  analyzed  in  periods  of  1  hr  maximum.  For  example, 
if  the  duration  were  4  hr,  it  would  be  treated  statistically  as  four  1-hr  cases.  If 
the  duration  were  4  hr  45  min,  it  would  be  treated  as  four  events  of  1  hr  and  one 
event  of  45  min.  (The  >  I -hr  time  periods  are  discussed  separately  later  in  the 
analysis.  ) 

When  the  measurements  from  the  31  days  were  subjected  to  these  limitations, 

185  events  with  time  periods  between  5  and  60  min  were  designated  as  the  data 
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base  for  analysis.  Each  event  consisted  of  a  30-sec  averaged  rain-rate  reading 
for  each  3 -sec  interval  throughout  the  event  period.  The  number  distribution  of 
the  185  events  with  respect  to  time,  mean  rain  rate  (R),  and  maximum  rate  IR^) 
are  shown  in  the  histograms  of  Figures  1  through  3. 


3.  DATA  ANALYSIS 


The  objective  of  the  investigation  was  to  use  these  high-resolution,  rain-rate 
data  in  an  attempt  to  discover  and/or  devise  the  means  of  predicting  the  variabil¬ 
ity  of  rainfall  from  standard  weather  observations  or  from  climatological  records. 

Variability  can  be  loosely  classified  by  describing  the  rain  type  as  widespread 
or  homogeneous  with  little  variability,  showery  with  moderate  variability,  and 
rates  that  vary  considerably  as  in  thunderstorms.  A  more  precise  measure  is 
that  of  standard  deviation  (a),  a  value  that  is  calculated  from  the  3-sec  rates  for 
each  particular  event.  When  o  is  multiplied  by  100  and  divided  by  the  mean  rate, 
it  becomes  a  relative  measure  where  a  is  expressed  as  the  percentage  of  the  mean 
value  as 
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Figure  4.  Number  Distribution  of  the  Coefficients 
of  Variation 


considerably.  It  thus  became  apparent  that  any  rate -prediction  scheme  would 
have  to  be  expressed  as  the  relationship  (Rr)  of  the  rates  (R)  to  their  means  as 


with  the  maximum  relative  rate  being 
R 

R  =  -21  . 

rm 

An  estimate  of  the  Rrm  value  for  a  given  V  can  be  made  by  using  the  linear 
regression  relationship  of 

«/R  =1  +  0. 0109  V  . 

rm 

The  plot  in  Figure  5  shows  the  v  Rrm  and  V  values  from  the  185  events.  The 
dashed  lines  indicate  the  standard  error  of  estimate  (SEE)  limits  of  ±  0.074.  If 
Eq.  (4)  is  solved  for  V  as 


s  *, 


iRrm-  I  +0  0109  V 


- SEE  =  0.07 
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Figure  5.  Plot  of  V  Rrm  vs  the  Coefficients  of 
Variation  for  the  185  Events 


we  now  have  the  means  to  determine  V  from  knowledge  of  R^  and  R. 

Rr  and  Rrm  as  defined  in  Eqs.  (2)  and  (3)  are  dimensionless  ratios  and  are 
used  as  such  in  the  equations  throughout  this  report.  However,  they  can  be  con¬ 
sidered  as  factors  of  the  mean  value  as  R  =  R  R  and  R  =  R  R  with  R  having 
units  of  mm  hr  1  when  comparing  relative  and  absolute  rates  as  used  in  the  fol¬ 
lowing  explanation  of  the  analytical  concept.  (Rain  rates  less  than  their  mean 
values  are,  in  all  probability,  not  germane  to  attenuation  problems  and  only 
Rr's>  1  are  considered  in  these  analyses.) 

The  coefficient  of  variation,  in  that  it  is  the  relative  measure  of  o,  is  a  rig¬ 
idly  defined  mathematical  quantity.  When  applied  to  a  normal  bell-shaped  curve, 
68  percent  of  the  points  in  the  data  set  will  be  included  in  plus  and  minus  one 
standard  deviation  from  the  mean  value.  Thus,  if  we  consider  a  theoretical  rate 
distribution  with  a  V  =  50  that  forms  a  normal  curve  with  R  =  1.  5  R  being  the  up¬ 
per  a  limit,  the  relative  rate  can  be  expressed  as  Rr  =  (R  +  a)/ R  or  1.  5.  In  this 
distribution,  the  values  of  R  =  0  to  1.5  R  contain  84  percent  of  the  points,  where¬ 
as  the  points  >  R  =  1.  5  R  are  16  percent  of  the  total  number.  Any  value  of  R  can 
similarly  be  defined  up  to  Rrm  and  the  percentage  larger  than  the  level  can  be 
calculated.  Figure  6  illustrates  this  principle  and  shows  several  Ry  levels  and 
the  percentages  of  points  that  are  above  those  values. 

Experimental  data  seldom  form  normal  distributions  thus,  some  deviation 
from  this  rule  can  be  expected.  However,  if  one  compares  rain-rate  events  hav¬ 
ing  the  same  V,  one  would  expect  reasonable  agreement  in  the  percentage  of  points 
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RAIN  RATE  (mm  hr'1 ) 


Figure  6.  Theoretical  Curve  Showing  Number  Distribution 
with  Rain  Rates.  The  dashed  vertical  lines  indicate  ±  one 
standard  deviation  and  the  Rr  relative  relationships  are  in¬ 
dicated  by  the  solid  line 

within  the  ±  a  limits.  Continuing  this  line  of  reasoning,  one  should  also  find  a 
comparable  percentage  of  points  above  any  level  that  is  defined  relative  to  the 
mean.  Since  each  point  of  a  rain-rate  plot  from  our  instrument  represents  a 
3 -sec  time  interval,  it  follows  that  the  percentage  of  time  that  the  rates  are  2 
any  specific  Rr  level  should  also  be  comparable. 

Figure  7  shows  1  hr  of  rain  rates  that  were  determined  from  weight  measure 
ments  obtained  at  HAFB  on  13  August  1982.  (Because  of  the  30-sec  averaging, 
the  total  time  is  actually  3570  sec  or  0.  992  hr. )  The  calculated  V  for  this  event 
is  50. 1.  The  horizontal  lines  represent  various  levels  of  Rr,  where  a  value  of  1 
is  the  mean  of  the  data.  A  trace  along  any  of  these  lines  will  delineate  the  times 
when  the  rate  becomes  equal  to  and  rises  above  that  level  and  another  when  it 
falls  below.  A  complete  line  scan  results  in  one  or  more  time  periods  when 
R  2  Rr  and  the  summation  of  the  periods  give  the  total  amount  of  time.  When 
time  is  expressed  as  a  percentage  of  the  total  time  of  the  event,  it  allows  compar 
isons  to  be  made  between  events. 

Table  1  lists  the  specific  time  periods  that  the  rates  in  Figure  7  were  2  the 

various  R  levels.  The  differences  in  these  values  and  those  from  the  normal 
r 

curve  of  Figure  6  are  the  effects  of  the  "non-normal",  natural  distribution  as 
presented  in  Figure  8. 


2130  2140  2150  2200  2210  2220  2230 

TIME 

Figure  7.  One  Hour  of  Rain  Rates  Having  a  Coefficient  of  Varia 
tion  of  50. 1.  The  horizontal  lines  indicate  various  levels  of  Rr 


RAIN  RATE  (mm  hr'1  ) 


Figure  8.  Number  Distribution  of  Rain  Rates  From 
the  Data  Plotted  in  Figure  7 


Table  1.  Time  Periods  Where  Rain  Rates  >  Various  Rr 
Levels  for  1  hr  of  Continuous  Rain  on  13  August  82 


13  August  82 

R  =  1.97  Rm  =  4.79 

Time  2130  -  2230  (3670  sec) 

Rrm  =2.44  V  =  50.  1% 

Time 

sec 

%  of  Total 

2138.3  -  2139.2 

50.6 

1.4 

2139. 6  -  2153. 0 

800.8 

22.4 

2155.2  -  2201. 5 

387. 8 

10.  9 

2206. 8  -  2211. 8 

1537.  1  sec  or  43.  1%  >  ' 

297.  9 

3r  of  L  0 

8.3 

' 

2140.  6  -  2152.2 

691.  3 

19.4 

2156.2  -  2159. 3 

182.6 

5.  1 

2207. 5  -  2211.2 

1095.9  sec  or  30.7%  5:  ' 

222.0 

of  1.2 

6.2 

2141. 5  -  2148. 4 

413.  1 

11.6 

2148.  6  -  2150.  1 

92.7 

2.  6 

2150. 3  -  2151.2 

50.6 

1.4 

2156.4  -  2157. 8 

81.  5 

2.  3 

2158.  4  -  2158.  9 

25.3 

0.7 

2207.7  -  2210. 0 

800.  8  sec  or  22 . 4%  a  R 

137.7 

r  of  1.4 

3.9 

2141.  9  -  2148.  0 

368.  1 

10.  3 

2148.  8  -  2149.  1 

16.  9 

0.5 

2149.  6  -  2150.  0 

22.  5 

0.6 

2150.  5  -  2150.  8 

16.  9 

0.  5 

2207. 8  -  2209. 1 
500.2  sec  or  14%  S  Rr  < 

75.  9 

>f  1.6 

2.  1 

2142. 1  -  2143.  1 

56.2 

1.6 

2144.4  -  2147. 8 

255.7  sec  or  7.2%  >  Rr 

199.  5 

of  1.8 

5.6 

2144.6  -  2147.5 

171. 4  sec  or  4.8%  S  Rr 

171.4 

of  2.0 

4.8 

Figure  9  shows  the  rain  rates  for  a  41-min  shower  on  25  August  1982  and  the 
distribution  of  the  rates  are  shown  in  Figure  10.  The  V  derived  from  these  rates, 
50,  is  nearly  identical  to  that  of  Figure  7  although  the  R  of  that  case  was  2.7 
times  larger  and  the  rate  distribution  obviously  different.  However,  the  percent¬ 
ages  of  time  that  the  rates  were  2:  the  various  Rr  levels  as  listed  in  Table  2  are 
comparable  to  those  of  Table  1. 


RAIN  RATE  (mm  hr'1 


Figure  9.  Rain  Rates  From  a  41  Minute  Shower  Having 
a  Coefficient  of  Variation  of  50.  The  horizontal  lines  in¬ 
dicate  various  levels  of  Rr 


Figure  10.  Number  Distribution  of  Rain  Rates  From 
the  Data  Plotted  in  Figure  9 


Table  2.  Percentages  of  Time  That  Rain  Rates  Were  2: 
Various  Rr  Levels  for  a  41 -min  Shower  on  25  August  82 


25  August  82  Time  0857:30 

Rr  =  0.73  Rm  =  1.58  Rrrn 

-  0938:30  (2456  sec) 

=  2.  16  V  =  50% 

Rr 

Time  (sec)  a  Rr 

%  of  Total  Time 

1.  0 

1109. 9 

45.2 

1.2 

716.  5 

29.2 

1.4 

539.  5 

22.  0 

1.6 

382.2 

15.6 

1.  8 

247.3 

10.  1 

2.0 

89.9 

3.7 

The  185  events  were  processed  in  the  identical  manner  as  described  for  those 
of  Figures  7  and  9.  The  percentages  of  time  that  the  rates  were  Rf  were  then 
plotted  vs  the  corresponding  V's.  Figure  11  shows  the  plots  for  Rr  =  1,  1.6,  and 
2.  5  The  solid  lines  indicate  the  trend  of  the  data. 

The  line  of  "best  fit"  for  Rr  =  1  over  the  range  of  V's  from  ~  20  to  ~160  can 
be  described  by 


T  =  52.7  e'°- 00474  V% 
P 


(6) 


where  Tp  is  the  percentage  of  time  that  the  rates  were  equal  to  or  larger  than  the 
mean  values.  If  the  rate  distributions  of  the  185  events  all  formed  normal  curves, 
T  would  equal  50  percent  for  each  case.  The  deviation  from  50  percent  can  be 
attributed  to  "non-normal"  or  skewed  distributions  (Figures  8  and  10)  and  the 
slope  of  the  line  marking  the  trend  in  the  upper  diagram  of  Figure  11  indicates 
that  the  skewness  increases  with  larger  V's.  The  trend  lines  in  the  other  two 
plots  show  a  changing  curvature  with  increasing  Rr- 

The  relationship  of  Tp  with  Rr  was  then  investigf ’.ed  and  it  was  found  that 


T  =  C  e*E  Rr  % 
P 


(7) 


for  ~  90  percent  of  the  R  values  or,  in  other  words,  for  R„  =  1  to  0.  9  R _ .  The 

r  r  r  rm 

remaining  ~  10  percent  of  Rr's,  from  Rr  =  0.  9  Rrm  to  Rrm»  have  percentages 

that  decrease  linearly  to  ~  0  at  Rrm-  These  relationships  are  shown  in  Figure  12 

for  the  rain  rates  of  Figure  7.  The  dashed  line  in  this  plot  is  the  semi-log  regres 

sion  line  derived  from  the  T  -  R  values  from  R  *  1  through  2.2.  The  solid  line 

p  r  r 

shows  the  linearly  interpolated  Tp's  from  Rr  *  2.2  to  Rrm  *  2.44. 


Figure  11.  Plots  Showing  the  Percentages  of  Time  That  the 
Rain  Rates  Were  >  Rr  of  1.  0,  1.  6,  and  2.  5  for  the  Coefficients 
of  Variation  of  the  185  Events.  The  data  trends  are  indicated 
by  the  solid  lines 


Although  scatter  was  evident  in  both  the  derived  E  and  C  values  of  Eq.  (7)  for 
the  ~  90  percent  Rr  grouping,  all  of  the  185  events  had  correlations  >  0.95  with  95 
of  them  >  0.  99.  When  E  was  plotted  vs  V  as  in  Figure  13,  a  log-log  relationship 
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- Tp  =522  .4  «  338  Rr 

- Tp  «  25.4  -  10.5  Rr 


X  T.  FROM  RAIN  RATES  OF  FIGURE  T 


Figure  12.  A  Plot  Showing  the  Tp  -  Rr  Mathematical  Rela¬ 
tionships  That  Describe  the  RainTlate  Data  of  Figure  7 


E  =  437.2  V 


described  the  data  in  the  best  manner  and  is  shown  by  the  solid  line.  The  dashed 
lines  represent  the  SEE  of  the  log  values,  ±  0.  096. 

The  values  of  C  from  the  185  events  for  the  ~  90  percent  Rr  grouping  are 
plotted  vs  V  in  Figure  14.  The  solid  line  in  this  plot  was  derived  by  substituting 
Eqs.  (6)  and  (8)  in  Eq.  (7)  for  Rr  =  1  and  solving  for  C  as 


C  =  52.7  e 


437.2  V'1-37  -  0.  00474  V 


The  SEE  (In)  for  these  data  was  0.  667 . 

When  Eqs.  (8)  and  (9)  are  used  in  place  of  E  and  C  in  Eq.  (7),  a  general  equa 
tion  incorporating  both  V  and  the  initial  90  percent  of  Rr's  can  be  derived  as 


T  =  52  7  437.2  V"1,37  (1  -  Rr)  -  .00474  V  (10 

P  e  %• 

The  interpolated  Tp  values  for  the  remaining  ~  10  percent  of  the  largest  Rr's 
are  described  by  the  equation 


Figure  13.  Exponents  From  the  Percentage  vs  Rr  Relation¬ 
ships  of  Eq.  (7)  for  Each  Coefficient  of  Variation  From  the 
183  Events.  The  siolid  line  was  derived  from  a  least-squared 
regression  analysis  and  the  dashed  lines  show  the  standard 
error  of  estimate 

437.2  V"1*37  (1  -  .9  R_)  -  .00474  V 

32.7  e  „ 

T 

P 


(Rrm  -  .9  Rr) 


’  70 


(1 


Figure  13  shows  estimates  T  vs  V  curves  for  various  R  levels  that  were 

p  r 

calculated  using  Eqs.  (10)  and  (II).  Figures  16  and  17  show  the  T^'s  that  were 

derived  from  the  rain  rates  of  Figures  7  and  9  compared  to  the  estimated  values 

calculated  using  Eqs.  (10)  and  (11).  [As  a  matter  of  interest,  if  Eq.  (10)  were 

applied  to  the  upper  a  limit  in  Figure  6  where  R  =  1.  5  for  V  =  50,  the  estimated 

T  =  14.9  percent  compared  to  the  actual  16  percent.  J  Table  3  lists  the  SEE  of 

tin  empirical  T^'s  from  the  185  events  for  several  Rr  levels  when  compared  to 

the  estimated  T  curves. 

P 

As  mentioned  previously,  there  were  16  days  with  continuous  rain  data  rang¬ 
ing  between  ~  1.  5  and  ~  5  hr.  Since  these  16  events  had  rain  rates  that  varied 
substantially  throughout  the  time  of  rainfall,  they  were  originally  divided  into 
periods  with  maximum  times  of  1  hr  as  explained  in  Section  2.  The  rates  from 
each  of  these  long-time  periods  were  now  processed  as  individual  events  to  test 


Figure  14.  Coefficients  From  the  Percentage  vs  Rr  Relation¬ 
ships  of  Eq.  (7)  for  Each  Coefficient  of  Variation.  The  solid 
line  was  derived  using  Eq.  (9) 

the  validity  of  the  equation  set  for  times  >  1  hr.  The  results  of  this  analysis  are 
compared  to  the  estimated  values  in  Table  4  and  show  that  the  statistics  that  de¬ 
scribe  the  1-hr  rain  events  are  also  applicable  to  times  >  1  hr. 

Although  Eqs  (10)  and  (11)  can  be  used  to  predict  the  amount  of  time  that  the 
rain  rates  exceed  specific  levels  of  intensity,  the  number  of  occurrences  during 
the  total  time  period  still  remains  a  question.  An  occurrence,  in  this  case,  is 
defined  as  a  period  of  time  starting  when  R  >  Rr  and  ending  at  R  <  Rr- 

When  one  considers  a  normal  curve  with  a  specified  V  as  in  Figure  6,  it  is 
apparent  that  the  total  number  of  points  in  the  distribution  can  be  raised  or  low¬ 
ered  (changing  the  values  of  the  y  axis)  with  no  effect  on  the  total  percentages  of 
time  for  the  various  Rr  levels  although  the  actual  number  of  points  will  be  differ¬ 
ent.  This  would  seem  to  indicate  that  the  number  of  occurrences  where  R  >  R 

r 

should  increase  in  accordance  with  longer  time  periods  for  events  of  the  same  V. 
However,  when  we  compared  six  events  of  similar  V's  (52  to  56)  that  ranges  in 
time  from  ~  9  min  to  —  5  hr,  we  found  vast  differences  in  the  number  of  occurren¬ 
ces  at  all  R  levels.  Table  5  lists  these  events  for  R„  =  1. 
r  r 

If  we  consider  Figure  6  once  again,  this  time  keeping  the  total  number  of 
points  constant  thereby  conserving  the  total  area  under  the  distribution  curve,  it 
is  again  apparent  that  the  profile  of  the  normal  curve  would  become  more  shallow 
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Figure  15.  Estimated  Percentages  of  Time  That  Rain  Rates  are 
Equal  to  or  Larger  Than  Various  Relative  Intensity  Levels  for 
Specified  Coefficients  of  Variation 


Table  3.  Standard  Error  of  Estimates 
From  the  Comparison  of  Tp's  From  the 
185  Events  to  the  Estimates  Calculated 
Using  Eqs.  (10)  and  (11) 


No.  Used  in 

R 

r 

SEE  % 

Analysis 

1. 0 

5.  56 

185 

1.2 

5.  42 

185 

1. 4 

4. 33 

183 

1.6 

3.  38 

174 

2.0 

2.  98 

137 

2.  5 

1.  98 

89 

3.  0 

1.  50 

56 

3.  5 

1.21 

40 

4.  0 

1.  15 

23 

5.  0 

0.  91 

10 

as  V  increases.  Thus,  one  might  deduce  that  for  identical  time  periods  the  num¬ 
ber  of  occurrences  at  any  specified  Rr  should  decrease  with  increasing  V's.  A 
comparison  of  events  with  identical  time  and  differing  V's  shows  a  slight  tendency 
towards  a  lower  number  of  occurrences  for  the  larger  V's  but  the  dispersion  is  so 
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Figure  16.  Comparison  Showing  Percentage  of  Time  That 
R  a:  Rr  From  Rain  Rates  of  Figure  7  and  Those  Calculated 
Using  Eqs.  (10)  and  (11) 


Figure  17.  Comparison  Showing  Percentage  of  Time  That 
R  S  Rr  From  Rain  Rates  of  Figure  9  and  Those  Calculated 
Using  Eqs.  (10)  and  (11) 


Table  4.  Standard  Error  of  Estimates 
From  the  Comparison  of  Tp's  From  16 
Events  >  1  hr  to  the  Estimates  Calcu¬ 
lated  Using  Eqs.  (10)  and  (11) 


Rr 

SEE  % 

No.  Used  in 
Analysis 

1.0 

4.  16 

16 

1.2 

4.43 

16 

1.4 

4.  07 

16 

1.6 

3.  55 

16 

2.0 

2.  98 

16 

2.5 

2.03 

14 

3.0 

1.33 

12 

3.  5 

0.  96 

8 

4.0 

0.74 

7 

5.  0 

0.29 

5 

Table  5.  Number  of  Occurrences  Where  R  >  Rr  -  1  From  Six  Events 
Having  Similar  Coefficients  of  Variation  and  Different  Times 


Date 

Total  Time 
hr  min  sec 

H 

No.  of 
Occurrences 

Tp  Where 

R  >  Rr  =  1 
% 

Tp  |Eq.  (10)] 

P  % 

9  August 

82 

9  19 

56.  1 

6 

44.  0 

40.4 

27  May 

83 

21  59 

52.  3 

15 

40.7 

41.  1 

12  October  83 

56  46 

56.  0 

5 

40.  8 

40.4 

13  August 

82 

1  48  31 

54.  1 

18 

38.  8 

40.  8 

20  July 

82 

3  20  29 

55.  1 

79 

39.  3 

40.  6 

16  July 

83 

4  58  28 

55.  9 

15 

46.  2 

40.4 

great  that  no  definite  conclusion  can  be  drawn.  Table  6  lists  12  events  of  1-hr 
duration  that  show  this  dispersion  at  the  Rf  =  1  level. 

Because  of  the  diversity  of  the  data  as  evidenced  in  Tables  5  and  6,  we  have 
been  unable  to  devise  a  means  to  reasonably  estimate  the  number  of  occurrences. 


4.  EFFECTS  OF  DATA  AVERAGING 


As  previously  mentioned,  the  rain  rates  used  in  the  preceding  analyses  were 
determined  from  3 -sec  weight  data  subjected  to  30-sec  averaging.  This  averag¬ 
ing  interval  was  found  to  be  optimum  for  our  particular  rain -rate  meter  in  that  it 


Table  6.  Number  of  Occurrences  Where  R  Rr  =  1  From  12  Events  of 
1-hr  Duration  With  Different  Coefficients  of  Variation 


Date 

H 

No.  of 
Occurrences 

Tp  Where 

R  >  Rr  =  1 

% 

Tp  [Eq.  ( 1 0) j 
% 

16  May  83 

27.5 

13 

46.6 

46.  3 

16  May  83 

31.  3 

6 

50.  5 

45.  4 

26  May  83 

32.4 

37 

48.  1 

45.2 

26  May  83 

41.  6 

41 

47.  0 

43.  3 

12  August  83 

45.  1 

5 

45.  1 

42.  9 

12  October  83 

50.  9 

1 1 

39.  7 

41.  1 

13  August  82 

53.  6 

4 

44.  1 

40.  9 

9  August  82 

63.  6 

10 

37.  8 

39.  0 

12  October  83 

83.  8 

5 

30.  7 

35.  4 

28  June  83 

96.  9 

5 

28.  9 

33.  3 

12  October  83 

127.  5 

5 

27.  8 

28.  8 

12  October  83 

133.  9 

1 

25.  7 

27.  9 

reduced  noise  while  retaining  the  characteristic  variability  of  the  rain  data. 
Averaged  rates  from  3 -sec  weight  data  are  found  by 


.  [ W(0  +  D)  yw(o  -h)J 
A  p  r\j 600 


(12) 


where  Av  is  the  weight  ±  n  time  increments  from  wn,  A  is  the  collection  area  in 

2  -3  u 

mm  ,  p  is  the  density  of  water  in  g  mm  ,  and  the  time  between  the  weight  read¬ 
ings  is  2  r)  X  3  sec/3600  sec  hr  1  or  n/600  hr.  When  wQ  advances  one  point  at  a 
time  through  the  3 -sec  weights  and  n  =  5,  a  30-sec  running-mean  (RM)  averaging 
results.  The  continuous  3 -sec  weights  may  also  be  used  to  average  in  the  more 
conventional  manner  when  Eq.  (12)  is  applied  as  w^  advances  in  increments  of  2r). 
Thus,  the  averaging  period  may  be  easily  changed  by  varying  rj- 

Figure  18  shows  the  30-sec  averaged  rates  as  in  Figure  7  and  the  determina¬ 
tions  from  a  10-min  RM  averaging  of  the  same  weight  data.  The  closed  circles 
are  the  six  readings  that  would  result  from  measurements  made  every  10  min.  As 
averaging  increases,  the  rate  variability  is  suppressed  and  substantial  differences 
may  be  present  in  the  values  of  the  minimum,  maximum,  and  V.  The  30-sec  RM 
values  divided  by  the  10-min  RM  in  this  figure  give  a  V  ratio  of  2.  36,  a  minimum 

rate  ratio  of  0.  44,  and  2.35  for  R  ratio. 

m 

The  determination  of  precipitation  rates  from  electronic -weight  measure¬ 
ments  has  a  decided  advantage  over  other  rate  measuring  methods  in  that  rate 
values  from  any  desired  averaging  interval  may  be  easily  calculated  from  the 
3 -sec  data  base.  In  comparison,  one  other  widely  used  instrument,  the  tipping- 
bucket  gauge,  records  when  the  bucket  fills  to  a  specified  capacity  and  dumps  the 


Figure  18.  Ten  Minute  Averaged  Rain  Rates  From  the  Weight  Data  Used 
in  the  30-sec  Averaged  Rates  of  Figure  7 


collected  water.  This  gives  time  periods  with  considerable  differences  in  light 
and  heavy  rain  situations.  Thus,  the  recorded  rate  data  from  these  type  measure¬ 
ments  are  already  averaged  over  varying,  non -consistent  time  intervals  with  long 
periods  in  light  rain  conditions  and  short-time  intervals  in  heavy. 

Accurate  averaging  may  be  of  considerable  importance  in  the  determination 
of  attenuation.  For  illustration,  we  can  consider  the  rate  profiles  in  Figure  18. 

The  rate  at  any  given  time  on  the  curves  in  this  figure  represent  the  amount 
of  rainfall  at  some  particular  point  on  the  earth's  surface  that  is  centered  along  a 
path  length  defined  by  the  wind  speed  and  averaging  period.  For  example,  if  the 
wind  velocity  during  this  hour  of  measurements  were  a  constant  1  m  sec  ,  any 
30-sec  rate  would  represent  the  rainfall  occurring  over  a  30-m  distance  in  the 
wind  direction.  The  10-min  averages  would  describe  600  m. 

We  may  assume  a  hypothetical  scenario  where  attenuation  measurements  were 
being  conducted  during  the  1-hr  rate  determinations  of  Figure  18.  We  may  further 
assume  that  the  path  length  was  a  known-fixed  distance  parallel  to  the  wind  direc¬ 
tion,  the  wind  speeds  were  being  recorded,  and  the  rates  were  obtained  from  a 
rain-rate  meter  (our  precipitation  weight-measuring  type)  situated  at  mid  range. 
With  these  stipulations,  it  is  apparent  that  the  weights  could  be  averaged  in  ac¬ 
cordance  with  the  wind  speed  to  produce  rates  that  correspond  to  path  length  dis¬ 
tances.  The  resulting  rate  structure  would  be  that  causing  the  attenuation. 


The  importance  of  averaging  caused  us  to  check  the  application  of  Eqs.  (10) 
and  (11)  that  were  derived  using  30-sec  data  for  the  estimation  of  Tp  to  those  rates 
determined  from  longer  averaging  periods. 

The  weight  data  that  produced  the  30-sec  rates  of  Figures  7  and  9  were  aver¬ 
aged  over  various  time  periods  from  1  to  30  min.  The  V's  and  R  's  resulting 
from  the  averaging  were  then  used  in  Eqs.  (10)  and  (11)  to  give  estimated  Tp's 
that  were  compared  to  the  measured  values.  The  comparisons  showed  reasonable 
agreement  for  averaging  intervals  up  to  10  min  as  evidenced  in  the  examples  of  the 
2 -min  rates  for  the  13  August  82  and  the  10  min  for  the  25  August  82  data  (Figure 
19).  The  Tp's  derived  from  >  10-min  averaging  did  not  conform  to  the  estimated 
curves  and,  at  this  time,  we  have  no  definitive  explanation  to  offer. 


Figure  19.  A  Comparison  of  Estimated  and  Measured  T„  Values  for  Two 
Averaging  Periods 


Thus,  for  rain  rates  averaged  £  10  min,  the  V,  R,  and  of  the  data  set 


may  be  used  for  estimating  Tp.  Rate  data  consisting  of  a  sufficient  number  of 


points  to  produce  a  valid  statistical  analysis  can  give  the  values  of  those  param¬ 
eters.  Hourly  or  6 -hr  weather  station  measurements  can  not  supply  this  informa¬ 
tion  nor  can  climatological  data.  But,  if  one  could  determine  V  from  knowledge 


of  averaging  time,  R  and  Rm  (Rrm),  values  that  may  possibly  be  deduced  from 


forecasts  or  climatology,  estimates  of  T  may  be  calculated.  We  thus  decided  to 


investigate  the  differences  in  the  signatures  of  rain  rates  averaged  over  various 
time  periods. 

Out  of  the  185  events  used  in  the  analysis,  63  were>  40-min  duration  and 
gave  a  sufficient  number  of  points  for  a  statistical  evaluation  from  30-min  RM 
averaging.  These  events  were  averaged  at  time  intervals  from  1  to  30-min  or 
from  n  =  10  to  300  in  Eq.  (12).  We  had  hoped  to  define  changes  in  V  and  Rrm  with 
averaging  time  or  with  respect  to  the  resulting  V  or  Rrm  from  averaging.  No 
correlation  was  found. 

However,  when  the  V's  and  Rrm's  were  compared,  they  formed  linear  rela¬ 
tionships  similar  to  the  30-sec  data  shown  in  Figure  5  as 

=  1  +  Ea  va  »  <13> 

where  the  subscript  "a"  refers  to  a  specific  averaging  period.  For  time  intervals 
^  10  min,  the  exponents  of  Eq.  (13)  (E&)  can  be  expressed  as 

Ea  =  0.0102  T*0- 089  min'1  (14) 

a  a 

with  a  SEE  (log)  =  0.  007  as  shown  in  Figure  20.  T&  is  the  time  of  the  averaging 
period  in  minutes.  This  plot  also  shows  the  non-conformity  in  the>  10-min  aver¬ 
aging. 

Substituting  Eq.  (14)  into  Eq.  (13)  yields  a  general  equation 

''  Rrma  =1  +  0.  0102  Ta'°‘  089  Va  (15) 

that  may  be  used  to  estimate  the  R  _  from  the  V  of  averaged  data  (s  10  min)  or  to 
*  rm 

find  the  estimated  V  of  the  data  set  when  R  is  known  by 

rm 

•J  RrmQ  -  1 

V  =  - V-089  %  •  (16) 

a  0.0102  Ta  °' 089 

When  Eq.  (16)  is  applied  to  the  10-min  RM  averaged  data  of  Figure  18  in 
which  Rrm  =  1.71,  the  estimated  V  =  36  vs  the  actual  37.  1.  Using  the  six  points 
from  the  conventional  10-min  averaging,  the  Rrm  =  3.  48/ 1.  96  =  1.78  and  the  es¬ 
timated  V  =  39.  1. 


VS; 


’.V 


E0»  0.0102  Ta  '0  089 


~+~ 


AVERAGE  TIME-  Ta  (min) 

Figure  20.  A  Plot  of  the  Exponents  of  Eq.  (13)  for  Various 
Averaging  Intervals 


5.  FINDINGS 

Intensities  of  rainfall  may  vary  considerably  from  one  rain  situation  to  an¬ 
other  giving  rain -rate  profiles  showing  small  changes  over  long-time  periods  or 
large  changes  occurring  within  a  few  seconds.  Since  many  endeavors  are  affected 
by  relatively  short  periods  of  rainfall  above  specific  levels  of  intensity,  such  as 
the  attenuation  experienced  by  electro -optical  sensing  systems,  the  ability  to  pre¬ 
dict  some  measure  of  the  rate  variability  from  forecast  information  could  prove 
essential.  Also,  evaluations  of  system  performances  in  given  geographical  areas 
may  be  possible  if  estimates  of  rate  intensities  could  be  deduced  from  statistical 
climatological  information. 

The  rain  investigation  detailed  in  this  report  was  conducted  with  the  objective 
of  establishing  a  method  to  predict  the  durations  and  frequency  of  specified  rate 
intensities  from  limited  inputs  such  as  those  obtained  from  forecasts  or  clima¬ 
tology.  Our  data  consisted  of  3  1  days  of  measured  rain  rates  from  which  185 
separate  events  were  isolated.  These  events  were  defined  and  classified  as  to 
their  rate  variabilities  in  terms  of  the  coefficient  of  variability  (V).  We  found 
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that  V  had  no  dependence  upon  the  amount  of  rainfall  that  occurred  but  did  change 
predictably  with  the  relative  measure  of  the  ratio  of  the  rates  (R)  over  the  mean 
values  (R)  as  R/R‘  Analysis  of  the  185  events  using  various  relative  rate-inten¬ 
sity  levels  showed  that  the  time  the  rates  were  equal  to  or  exceeded  the  relative 
levels  varied  with  V  and  could  be  mathematically  modeled  or  described  by  an 
equation  set.  We  also  found  that  the  relationships  of  the  relative  maximum  rates 
and  their  corresponding  V's  varied  predictably  with  averaging  period  up  to  10-min 
time  intervals. 

Thus,  we  have  deduced  a  method  by  which  rain  variability,  in  terms  of  per¬ 
centage  of  time  the  rates  are  >  intensity  levels  relative  to  their  mean  values,  may 
be  estimated  from  knowledge  or  deduction  of  two  parameters,  the  mean  and  max¬ 
imum  rates. 

For  example,  if  one  could  obtain  the  mean  and  maximum  rain  rates  that  have 

occurred  or  are  forecasted  to  occur  during  a  specified  time  period,  one  may  apply 

the  relative  maximum  rate  (R  )  and  time  interval  (T  )  in  Eq.  (16)  to  estimate  the 

rm  a 

coefficient  of  variation  for  that  particular  situation.  The  estimated  percentage  of 
total  time  that  the  rain  rates  are  equal  to  or  become  larger  than  any  relative  value 
(Rr),  where  Rr  is  the  factor  above  the  mean,  may  then  be  determined  from  the 
plots  of  Figure  15  or  calculated  using  Eqs.  (10)  and  (11).  In  this  scenario,  one 
might  find  it  more  convenient  to  estimate  Tp  using  the  isolines  of  V  from  the  plot 
of  Tp  vs  R^  from  Figure  21. 

On  the  negative  side,  we  failed  to  devise  the  means  to  estimate  the  frequency 
of  occurrence  for  rates  rising  above  relative  levels.  Neither  could  we  find  a 
method  to  relate  the  change  in  variability  from  different  averaging  periods. 

Although  the  equations  in  this  report  have  been  developed  for  averaging  per¬ 
iods  s  10  min,  the  data  from  longer  time  intervals  may  be  used  to  provide  less 
accurate  values  of  percentage  of  time  if  one  takes  into  account  that  the  estimated 
results  may  be  low  by  as  much  as  25  percent  at  the  relative  level  of  1  to  60  per¬ 
cent  when  approaching  the  relative  maximum  level. 

Finally,  we  recognize  the  fact  that  the  rain  data  used  in  the  preceding  analy¬ 
ses  were  obtained  at  one  geographical  location  and  do  not  necessarily  represent 
the  particular  characteristics  found  in  other  global  areas.  However,  since  we 
classified  each  rain  situation  using  a  mathematical  definition  of  relative  variability 
and  used  relative  values  of  rain  rates  in  the  analyses,  we  believe  that  the  equa¬ 
tions  presented  herein  may  be  applicable  in  other  rain  situations  regardless  of 
location. 
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